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ABSTRACT 



We present spectrophotometry in the 3600-9700 A region for a sample of 39 HII 
regions in the Galaxy and Magellanic Clouds, for which independent information is 
available on the spectral types and effective temperatures of the ionizing stars. The 
spectra have been used to evaluate nebular diagnostics of stellar temperature, metal 
abundance, and ionization parameter, and compare the observed behavior of the line 
indices with predictions of nebular photoionization models. We observe a strong degen- 
eracy between forbidden-line sequences produced by changes in stellar %ff and metal 
abundance, which severely complicates the application of many forbidden-line diagnos- 
tics to extragalactic HII regions. Our data confirm however that the Edmunds & Pagel 
[O II] + [0 III] abundance index and the Vflchez & Pagel rf index provide more robust 
diagnostics of metal abundance and stellar effective temperature, respectively. A com- 
parison of the fractional helium ionization of the HII regions with stellar temperature 
confirms the reliability of the spectral type vs %jj calibration for the relevant temper- 
ature range %ff < 38000 K. We use empirical relations between the nebular hardness 
indices and T^g to reinvestigate the case for systematic variations in the stellar effective 
temperatures and the upper IMFs of massive stars in extragalactic HII regions. The 
data are consistent with a significant softening of the ionizing spectra (consistent with 
cooler stellar temperatures) with increasing metal abundance, especially for Z < Zq. 
However unresolved degeneracies between Z and %g still complicate the interpretation 
of this result. 
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1. INTRODUCTION 

Over the past 15 years high-quality emission-line spectra have been obtained for hundreds of 
HII regions in nearby spiral and irregular galaxies. The primary use of these data has been to 
study the chemical composition patterns in galactic disks (e.g., McCall, Rybski, & Shields 1985, 
Vila-Costas &z Edmunds 1992, Zaritsky, Kennicutt, & Huchra 1994, Dinerstein 1996). However 
with the availability of state-of-the-art photoionization codes such as cloudy (Ferland etal. 1998), 
combined with stellar population synthesis codes for young clusters (e.g., Leitherer & Heckman 
1995, Leitherer et al. 1999), HII region spectra are being applied increasingly to measure the ages 
and initial mass functions (IMFs) of the ionizing star clusters (e.g., Vflchez & Pagel 1988, Shields 
1990, Stasihska & Leitherer 1996, Bresolin, Kennicutt, & Garnett 1999, hereafter BKG). Similar 
techniques have been applied to model the infrared emission-line spectra of IR-luminous starbursts, 
and the results have been interpreted as evidence for an anomalous IMF in these objects (e.g., Rieke 
et al. 1993). 

A basic limitation in this approach is its reliance on a long chain of theoretical inputs: stellar 
evolution models as functions of stellar mass, chemical composition and mass-loss rates; stellar 
atmosphere models as functions of effective temperature, surface gravity, chemical composition, 
and mass-loss properties, and photoionization models for the surrounding H II region. The derived 
nebular abundances tend to be insensitive to the details of these models, but the inferred properties 
of the ionizing stars are sensitive to the model inputs at every step in the chain. The problem is 
especially acute when using H II region or starburst spectra to constrain the effective temperatures 
and IMFs of the ionizing stars (e.g., Mathis 1985, McCall et al. 1985, Vflchez k Pagel 1988, BKG). 
Most studies show evidence for a systematic softening of the ionizing continuum of HII regions 
with increasing metal abundance, but the absolute range in stellar effective temperatures is highly 
model dependent, and the case for a systematic variation in the upper IMF is shaky at best (BKG 
and references therein). 

Another approach, which can circumvent many of these difficulties, is to test the nebular 
diagnostics empirically, by obtaining spectra for nearby H II regions which are ionized by stars of 
known spectral type. This makes it possible to anchor the nebular indicators of stellar temperature 
directly to the stellar spectral classification system. Furthermore, by modeling these HII regions 
with the same methods that are applied to extragalactic regions one can directly assess the reliability 
of the model-based approach. 

This empirical approach was first explored by Chopinet & Lortet-Zuckermann (1976) and Kaler 
(1978) to calibrate the [O IIl]A5007/H/3 ratio as an indicator of stellar effective temperature. For 
studies of extragalactic H II regions and starbursts a more robust hardness index than [O III] /H/3 
is desirable, because the excitation of [O III] varies locally within HII regions, and it is sensitive 
to other physical parameters (e.g., metal abundance, ionization parameter, dust) which can vary 
systematically in galaxies (e.g., Kennicutt 1984, Shields 1990). Other hardness indices have been 
developed for extragalactic applications (e.g., Vflchez & Pagel 1988), but their reliability has not 
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been tested empirically. 

In this paper we report the results of a spectrophotometric study of 39 HII regions in the 
Galaxy, the LMC, and the SMC, which is aimed at providing an empirical foundation for nebular- 
based measurements of the stellar ionizing continuum and IMF in galaxies. We use the spectra 
to evaluate several nebular diagnostics of massive stellar populations, and to compare the known 
properties of the stars in these regions with those inferred from photoionization modeling. Our 
study is mainly motivated by applications to stellar temperatures and IMFs, but the approach is 
more generally applicable to spectral modeling of extragalactic HII regions and starbursts. This 
paper is complementary in approach to a recent analysis of Oey et al. (2000), which has addressed 
many of these same questions using detailed point-by-point observations of 4 H II regions in the 
LMC. 

The remainder of this paper is organized as follows. The HII region sample and the obser- 
vations are described in § 2 and § 3, respectively. In § 4 we analyze the behavior of the principal 
diagnostic line ratios and investigate several commonly applied abundance and spectral hardness 
indicators. In § 5 we apply the empirical calibration to constrain the range of stellar effective 
temperature and upper IMF in galactic disks. 

2. H II REGION SAMPLE AND PROPERTIES 

Since the first prerequisite for this work is knowledge of the exciting stars in the calibrating 
H II regions, the ideal approach is to observe small Galactic H II regions ionized by single stars (or 
a few stars) of known spectral type. Such "calibrating" objects comprise approximately half of our 
sample, while the remaining objects, ionized by larger OB associations or clusters, provide the link 
to the brighter class of H II regions observed in external galaxies. 

The calibrating HII regions were selected to span the maximum available range in stellar 
spectral type and effective temperature. Preference was given to nebulae with angular diameters 
that fit within the length of the spectrograph slit (3' — 6'), and to regions with low to moderate 
extinction, with E(B — V) < 2 mag in most cases. Most of our sample was drawn from the northern 
survey of small Galactic H II regions defined by Hunter & Massey (1990; =HM90), and the southern 
sample of HII regions studied by Shaver et al. (1983) in their abundance survey of the Milky Way. 
Identifications and spectral types for the ionizing stars were taken directly from HM90 for the 
northern sample, and from a variety of published sources for the southern objects. Most of the 
latter derive from classification work in the 1970's by Georgelin and collaborators (e.g., Georgelin, 
Georgelin, & Roux 1973, Chopinet, Georgelin, & Lortet-Zuckermann 1973). Table 1 lists the main 
properties of the H II regions in our sample, including the sources of the stellar spectral types and 
oxygen abundances. Column 8 indicates the observatory used to obtain the spectra, and asterisks 
denote objects in the calibrating subsample. 

It is relatively easy to compile an adequate sample of small H II regions ionized by single stars 
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of type 06V and later (%ff < 42000). However this becomes nearly impossible for hotter O stars, 
because such stars almost always form in rich OB associations along with numerous later-type 
stars. The contribution of the cooler stars to the ionizing continuum can be significant, and this 
composite ionization needs to be taken into account in the determination of %g. We observed 
several of these large Galactic HII regions, including M8, M16, M17, the Rosette nebula (S275), 
the Carina nebula, and the giant H II region NGC 3603. We used the most recent studies of their 
stellar contents, including near-infrared observations in many cases, to identify all of the principal 
ionizing stars. In order to extend our coverage to even higher stellar temperatures, we also observed 
four H II regions ionized by early-type Wolf-Rayet stars, the Galactic regions RCW 5 = NGC 2359 
(WN4), RCW 48 = NGC 3199 (WN5), and the LMC regions DEM 174 (WN4p) and DEM 231 
(WN3). 

The remainder of our sample consists of much larger HII regions located in the LMC and 
SMC. Most of these objects are physically distinct from the smaller Galactic HII regions, being 
ionized by much more populous OB associations, and as such are less suitable for our empirical 
calibration of nebular diagnostics. However these objects are directly analagous to the more distant 
extragalactic H II regions, and their stellar contents have been cataloged by various authors (Table 
1), so we can analyze them on the same basis as the smaller Galactic HII regions. 

To assign a stellar effective temperature to the HII regions we have adopted the spectral 
type-2^ calibration by Vacca, Garmany and Shull (1996), together with the stellar ionizing fluxes 
predicted by Schaerer & de Koter (1997). The resulting temperatures and luminosities are listed 
in columns 4 and 5 of Table 1. In the case of an ensemble of hot stars, a mean temperature was 
derived by summing the He and H° ionizing fluxes of all of the stars, and calculating the effective 
temperature of a class V star having the same Qi/Qo ratio (He to H° total ionizing flux ratio). 
Effective temperatures for the W-R stars were estimated using the semi-empirical calibration of 
Esteban et al. (1993). More detailed information on the stellar contents of the individual HII 
regions is given in the Appendix. 

Oxygen abundances, expressed as 12+log(0/H), are given in column 6, and their sources are 
indicated in column 7. When abundances were unavailable in the literature and could not be 
measured directly from our spectra, we estimated the abundance using the empirical indicator 
R 23 =([0 III]AA4959,5007+[O II]AA3726,3729)/H/3 (Pagel etal. 1979), as calibrated by Edmunds 
and Pagel (1984). These estimates are listed in parentheses in Table 1. The accuracy of the 
i?23-based abundances is of order ±0.2 dex, but this is sufficient for our purposes. 

Figure 1 shows the range of stellar temperatures and oxygen abundances for the H II regions in 
our sample. Squares and round points denote calibrating and non-calibrating HII regions with di- 
rect (electron temperature based) abundance measurements, respectively. Open triangles indicate 
those regions with empirical (R23) abundance estimates. The HII regions cover a stellar tempera- 
ture range of 32000 - 48000 K, with the Wolf-Rayet nebulae extending to 67000 K. The 28 Galactic 
HII regions span a range of O/H abundances of —0.5 < [O/H] < 0.0 (assuming 12+log(O/H) = 
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8.9), while the 9 LMC regions cluster around [O/H] ~ —0.4, and the single SMC region (N66 = 
NGC 346) has [O/H] = -0.7. 

The range of abundances in the sample is important for the interpretation which follows, 
because changes in metal abundance can mimic many of the trends in nebular spectra that are 
produced by changes in stellar temperature or IMF. For lack of a better prescription we have 
adopted a single conversion of spectral type to %g in our analysis, independent of abundance, and 
this conceivably could introduce some coupling between the two parameters (§ 4.3). Fortunately, 
Fig. 1 shows that Z/Zq and %g are largely decoupled from each other in our sample, and this 
should allow us to separate their effects on the nebular spectra. 

3. OBSERVATIONS AND DATA REDUCTION 

Our data consist of spectrophotometry covering the wavelength range 3600-9700 A. The ex- 
tended coverage to the near-infrared includes the [S III] AA9069, 9532 doublet, which when combined 
with measurements of [S II] AA6716, 6731, [O II] A3727, and [O III] AA4959, 5007 provides a robust 
measure of the hardness of the ionizing stellar continuum (Mathis 1985, Vflchez & Pagel 1988). 

3.1. Steward Observatory Data 

Spectra of the northern sample of Galactic H II regions were obtained in 1993 Oct and 1995 Jun 
with the B&C CCD spectrograph on the Bok telescope, equipped with a thinned Loral 800x1200 
element CCD detector. Two grating settings were used to cover the full spectral range. A 400 
g mm -1 grating blazed at 4800 A provided coverage of the 3600-6900 A region, hereafter denoted 
as the "blue" spectra. A second set of "red" spectra were obtained with a 400 g mm grating 
blazed at 7500 A to cover the 6400-9700 A region. A slit width of 2 / .'5 provided a spectral resolution 
of 8 A FWHM. The full spatial coverage of a spectrum was 3f5, with a single pixel on the detector 
projecting to / /83 on the sky. Spectra of NGC 7538 and NGC 7635 were obtained in in 1988 Nov 
using the same spectrograph but with a 800x800 element Texas Instruments CCD detector. The 
same setup was used for the red spectra, but the blue spectra were obtained with a 600 g mm 
grating blazed at 3570 A, to cover the region between 3700-5100 A. 

Figure 2 shows images of the HII regions, with the slit positions superimposed. Most of the 
images were taken from the Digitized Sky Surveys^. The HII regions were observed at a fixed 
position extending E-W through the center of the nebula (but offset from bright stars to avoid 
contamination of the spectrum), and also in a drift scanned mode, where the image of the HII 
region was trailed across the slit to provide an integrated spectrum. Comparison of the fixed 
and scanned spectra (integrated over the length of the slit in both instances) reveals surprisingly 
small differences in the main diagnostic line ratios, indicating that the fixed pointings provide a 
representative sampling of the nebula. This is illustrated in Figure 3, which compares the behavior 
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of [O IIl]/H/3 vs [N II] /Ha and [S II] /Ha for the two sets of spectra, with each pair of observations 
connected. Most of the differences follow the trajectories expected from changes in the nebular 
ionization parameter. We have chosen to analyze the fixed position data, in order to take advantage 
of their higher signal/noise, and to maintain consistency with the CTIO data, which were taken 
as fixed pointings. The only exception is the very large HII region M8 (S25), which shows a large 
inconsistency between fixed and scanned spectra, presumably because the former is is strongly 
influenced by local variations in ionization structure. In that case we use the drift scanned spectrum, 
which samples a region of 10' x 3' (EW x NS), centered on the position listed in Table 1 (see 
Fig. 2d). 

Exposure times of 300—1200 sec were chosen to provide accurate spectrophotometry of the 
bright diagnostic lines (S/N > 20 in most cases). The fainter auroral lines such as [O III] A4363 
were not a target of this study, and in many cases they were too weak in our spectra to be useful. 
The spectrophotometric flux calibration was determined using standard stars from Massey et al. 
(1988) for the blue spectra, and stars from Oke & Gunn (1983) for the red spectra. The observations 
were made with 2"5 and 4"5 wide slits, oriented to limit errors from atmospheric dispersion effects 
to <0.05 mag across the wavelength range. 

3.2. CTIO Data 

The southern Galactic HII regions and all of the LMC and SMC HII regions were observed 
in 1987 Jan, 1987 Dec, and 1988 Jan using the CTIO 1.0 m and 1.5 m telescopes. At the time 
that these observations were made a blue-sensitive CCD detector was not available, so the blue 
spectra were obtained with a 2D FRUTTI photon counting detector system on the 1.0 m Cassegrain 
spectrograph. The instrument was configured with a 600 g mm -1 grating blazed at 5000 A and 
WG-360 blocking filter, which provided coverage of 3600-7000 A and a resolution of 6 A FWHM 
when used with an 8" x 6' slit. All observations were made with the slit oriented E-W at a 
fixed position, usually centered 10" south of the central star. Due to the limited dynamic range 
of the 2D FRUTTI detector, observations of bright regions were made with neutral density filters, 
and exposure times for most objects were 900-1800 sec, largely independent of surface brightness. 
Standard stars from Stone & Baldwin (1983) were observed with a 12" wide slit for flux calibration. 

Red spectra were obtained with the Cassegrain spectrograph on the CTIO 1.5 m telescope, 
equipped with an unthinned 385x576 GEC CCD detector. Two grating settings were used for most 
of the objects. A 158 g mm -1 grating blazed at 8000 A in first order and a OG-550 blocking filter 
provided complete coverage of the 5730-10200 A region, with a resolution of approximately 20 A 



8 The Digitized Sky Surveys were produced at the Space Telescope Science Institute under U.S. Government grant 
NAG W-2166. The images of these surveys were based on photographic data obtained using the Oschin Schmidt 
Telescope at Palomar Mountain and the UK Schmidt Telescope. The plates were processed into the present digital 
form with the permission of these institutions. 
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FWHM when used with a 6" x 6' slit. A second set of spectra were obtained using a 600 g mm 1 
grating blazed at 6750 A, to cover the region 5650-6850 A at 5 A resolution, again with a 6" slit. 
The low-resolution spectra were used primarily to measure the diagnostic lines longward of 7000 
A (mainly [S III] AA9069, 9532 and [Ar III] A7135), while the higher-resolution spectra provided the 
best measurements of the [N II] AA6548, 6583, [S II] AA6716, 6731, and He I A6678 lines. Exposure 
times ranged from 120 to 1800 sec. Standard stars from Stone and Baldwin (1983) and Baldwin & 
Stone (1984) were observed with a 12" slit. 

In order to investigate the local variation in the nebular diagnostic indices within individual 
HII regions, we obtained spectra at additional positions in the Carina and 30 Doradus HII regions, 
as shown in Figures 2c and 2g, respectively. The observations of 30 Dor include a series of five 
positions extending E-W a total of 24' through the center of the HII region. We also obtained 
red spectra only at 22 other positions in 30 Dor and 12 positions in the Orion nebula, to match 
previous observations at 3700-7200 A by Mathis, Chu, & Peterson (1985) and Peimbert & Torres- 
Peimbert (1977), respectively. Charts showing the positions of these measurements can be found in 
the original papers. The red spectra were used to extend the wavelength coverage of the published 
measurements to the [S III] AA9069, 9532 lines. We used the overlapping coverage of the Ha, [N II], 
and [S II] lines to firmly tie in the [S III] fluxes, and confirm that the identical regions were observed. 



3.3. Data Reduction and Calibration 

The spectra were reduced using the two-dimensional spectrum reduction package in IRAFQ, 
following conventional procedures. We limit the discussion here to reductions and tests that were 
unique to this data set. Spatial distortions in the long-slit spectra were removed using comparison 
lamp exposures. The distortions were considerable in the 2D FRUTTI data, and data within 5-10% 
of the slit edges were excluded in subsequent reductions. The 2D FRUTTI spectra were checked 
for deadtime effects in the bright lines, using short exposures obtained with darker neutral filters, 
and by comparing the observed [O III] AA4959, 5007 doublet ratio (the strongest feature in most 
raw spectra) to its theoretical value. Deadtime effects were present in only a small fraction of the 
spectra, and in such cases the [O III] A5007 line was discarded, and the scaled flux of [O III] A4959 
was used instead. Similar checks were carried out with multiple neutral density filters on the 
standard stars, and measurements with any hint of a saturation problem were discarded. 

Calibration of the red CCD spectra in the [S III] AA9069, 9532 region also required special 
care, due to the presence of unresolved telluric water absorption lines, To maximize the reliability 
of the [S III] measurements we made frequent observations of subdwarf standard stars located at 
the same airmass as the HII regions, and used a high-order spline fit to the stellar spectra to 



9 IRAF is distributed by the National Optical Astronomical Observatories, which are operated by AURA Inc., 
under contract with the National Science Foundation. 
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remove the absorption on wavelength scales of > 50 A (exclusive of the regions around the Paschen 
emission lines, which are influenced by absorption lines in the standard stars). The reliability of this 
procedure was confirmed by the degree of consistency of the observed [S III] AA9069, 9532 doublet 
ratios with its theoretical value (2.44), and by the shapes of the calibrated continuum spectra in 
the HII regions themselves. As a conservative measure we attach a mean uncertainty of ±20% to 
the summed [S III] AA9069, 9532 fluxes. 

The overall precision of the spectrophotometry was determined from comparisons of multiple 
observations of the same objects, intercomparison of the blue and red spectra in the overlapping 
wavelength regions, and from the measured ratios of the [O III] and [S III] doublets. The typi- 
cal uncertainties for strong lines range from ±5% or better for the Bok CCD spectrophotometry 
(excluding the [S III] lines) to ±10% for the CTIO blue spectra (2D FRUTTI) and ±20% for the 
[S III] fluxes and the CTIO [O II] fluxes. External comparisons with published sources such as 
Dufour (1975), Shaver et al. (1983), Mathis et al. (1985), and Hunter (1992) show consistency 
to the degree that is expected given the measuring errors above and differences in positions and 
aperture sizes. The only exception is a tendency for the He I A5876 fluxes measured at CTIO (2D 
FRUTTI) to be significantly lower than those from published sources, and consequently we attach 
more weight to the He I A6678 (CCD) data in these objects. 



3.4. Emission-Line Measurements 

A one-dimensional spectrum of each region was extracted by summing over the radial extent of 
the HII region on the slit. Whenever possible a background sky measured from the outer parts of 
the slit was subtracted, but when the H II region filled the slit the sky was included and subtracted 
in the fitting of each line flux (after confirming the absence of sky line contamination). Analysis 
of the extracted ID spectra followed standard procedures. Line fluxes were measured by a direct 
integration of the line profile, with subtraction of a linearly interpolated continuum, using either 
the SPLOT task in IRAF or our own software. The flux scales of the blue and red spectra were 
tied together using the sum of the Ha and [N II] AA6548, 6583 lines. 

Reddening corrections were derived from the Balmer decrement (Ha, H/3, Hy), the theoretical 
Balmer line ratios as calculated by Hummer & Storey (1987), and the average interstellar reddening 
curve from Cardelli, Clayton, & Mathis (1989). When measured electron temperatures and densities 
were available, either from the literature or from our spectra, they were used to determine the 
theoretical Balmer decrement; otherwise a temperature of 10000 K and density of 100 cm -1 were 
assumed. 

Table 2 lists the reddening-corrected line fluxes for the principal diagnostic lines, as well as 
the logarithmic extinction at H/3 (column 2). All fluxes are normalized to f{H(3) = 100. For HII 
regions with multiple measurements we have listed the fluxes for a single representative position. 
Complete line lists, including the fainter lines and the multiple position measurements are available 
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from the authors upon request. 

4. ANALYSIS 

Before we specifically address the nebular diagnostics of abundance and stellar temperature, 
it is instructive to examine the behavior of the H II region spectra using the well-known diagnostic 
diagrams of Baldwin, Phillips, & Terlevich (1981), and compare them to the spectra of more distant 
and luminous extragalactic HII regions. 

To aid in the physical interpretation of these diagnostic diagrams, we have computed a set 
of cloudy nebular models (Ferland et al. 1998), using non-LTE stellar atmosphere models by the 
Munich group (Pauldrach et al. 1998) as the source of ionizing photons, as described in BKG. These 
will be used later to compare the stellar effective temperatures from nebular model fitting to the 
actual effective temperatures of the exciting stars. However in this section the main use of the 
models will be to illustrate how changes in metal abundance, stellar temperature, or ionization 
parameter affect the various line ratios. 

4.1. Diagnostic Diagrams and Ionization Parameter 

Figure 4 shows two of the most commonly used diagnostic diagrams for extragalactic HII 
regions and AGNs, with the excitation ratio [O IIl]A5007/H/3 plotted against [N II]A6583/Ha and 
[S II] AA6716, 6731/Ha. Here and in the following figures we have subdivided the sample: solid 
squares mark the calibrating Galactic HII regions ionized by one (or a few) OB stars; open squares 
correspond to the more luminous Galactic H II regions ionized by associations of OB stars with a 
range of spectral types; open circles denote the LMC and SMC regions, and open stars indicate the 
four nebular ionized by early-type W-R stars. For comparison we have plotted with small triangles 
the same line ratios for a large sample of extragalactic HII regions in Sa-Im galaxies, from BKG. 

Extragalactic HII regions show a tight spectral sequence in both of these diagrams (McCall, 
Rybski, & Shields 1985, Kennicutt & Garnett 1996), which often is interpreted as evidence that gi- 
ant HII regions are predominantly radiation bounded objects (McCall et al. 1985). Figure 4 shows 
that the much smaller Galactic HII regions also lie on this same sequence. The correspondence 
is especially tight in the case of [O IIl]/H/3 vs [N II] /Ha, where the sequences are virtually indis- 
tinguishable. The WR nebulae lie well above the sequence, which reflects the very high effective 
temperatures of these stars (T^jj > 55000 K). The LMC HII regions show a slight tendency toward 
weaker [N II] emission, which probably reflects the anomalous N/O ratio in this galaxy (Garnett 
1999). 

In extragalactic HII regions the tight spectral sequences (Fig. 4) are mainly interpreted as 
being abundance sequences (e.g., Searle 1971, Shields & Searle 1978, Pagel & Edmunds 1978, Mc- 
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Call et al. 1985). However in our sample the primary variable along the sequence is not abundance 
but stellar temperature. This is shown in Figure 5, which again shows the [O IIl]/H/3 vs [N II] /Ha 
relation, but this time with the points coded by stellar temperature (left) and by oxygen abundance 
(right). Although both parameters clearly influence the position of HII regions along the sequence, 
variations in %g clearly dominate in this sample. The virtual indistinguishability between this %g 
sequence with the extragalactic "abundance" sequence illustrates the degeneracy between temper- 
ature and abundance (and to some degree ionization parameter) in this type of diagnostic diagram. 
In hindsight this may help to explain why it has proven so difficult to disentagle variations in 
ionization temperature and IMF from abundance variations in extragalactic H II regions. 

The correspondence between the Galactic and extragalactic sequences is much poorer for 
[O IIl]/H/3 vs [S II] /Ha, with the small Galactic HII regions showing systematically weaker [S II] 
emission (Fig. 4). A comparison with the driftscanned spectra in Fig. 3 shows that this is partly 
due to the spatial undersampling of our slit spectra, but the remainder of the difference appears to 
be a systematic difference in mean ionization parameter between the two samples. This is quan- 
tified in Figure 6, which plots the ionization-sensitive ratio [S II] AA6716, 6731 /[S III] AA9069, 9532 
as a function of T^g of the ionizing stars (same symbols as Fig. 4). The The [S II] /[S III] ratio is 
commonly used as a diagnostic of the nebular ionization parameter U (Diaz et al. 1991), defined 
in this case as the ratio of ionizing photon to electron densities. It can be expressed as: 

U = -%* (1) 
AttRs nc 

where Q H o is the ionizing luminosity of the stars (photons per sec), R$ is the Stromgren radius of 
the HII region, n is the electron density, and c is the speed of light. Superimposed on Fig. 6 are 
photoionization models for values of \ogU between —1.5 and —4.0 (see BKG for details). Most of 
the regions studied here lie in the range of —3.5 < log U < —1.5, whereas most bright extragalactic 
HII regions lie in a narrower range with —3.5 < \ogU < —2.5 (BKG). The difference primarily 
reflects the larger aperture sizes in extragalactic studies (typically hundreds of parsecs in projected 
diameter). This difference is not large enough to affect the comparison of [/-insensitive indices 
such as the abundance parameter R23 and the nebular hardness parameter rj (§4.3), but it may 
need to be taken into account when applying more ionization-sensitive indices such as [O IIl]/H/3, 
[Ne IH]/H/3, or [O III]/[N II] (below). 



4.2. Abundance Indices 

The backbone of the extragalactic abundance scale rests on direct electron temperature based 
measurements, but in metal-rich H II regions the temperature-sensitive auroral lines are unobserv- 
able, and most abundance information in that regime is based on "empirical" excitation indices 
such as R 23 (Pagel et al. 1979) or [O III]A5007/[N II]A6583 (Alloin et al. 1979). Our data are not 
of sufficient quality to improve on the existing calibrations of these indices, but we can use our 
spatially resolved observations of 30 Dor, Carina, and Orion to investigate the robustness of the 
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empirical abundance indices. Do spectra of different regions in the same object yield consistent 
abundances? 

As a test of the R23 index, the left panel of Figure 7 shows the relation between [O III] AA4959, 5007/H/3 
and [O II] A3727/H/3 (with linear scales) for the multiple positions in 30 Dor (open stars), Carina 
(triangles), and Orion (open circles). The 30 Dor data include spectra from Mathis et al. (1985) 
and the Orion data are all from Peimbert & Torres- Peimbert (1977). All three HII regions show 
a very large range of excitation across the regions sampled, with both [O III]/H/? and [O Il]/H/3 
varying by factors of 2-5. However the sum of the line fluxes (-R23) is relatively constant, as shown 
by the diagonal lines in Fig. 7, The constancy of R23 is especially impressive in 30 Doradus and 
Carina, where the dispersion in R23 transforms to a full range of <0.1 dex in O/H. The R23 index 
is less well-behaved in Orion, though the range in abundances about the mean value is no larger 
than the calibration uncertainty of the method. Our results confirm a similar test applied to the 
giant HII region NGC 604 in M33 (Diaz et al. 1987), and suggest that even spatially undersampled 
observations of local HII regions may be useful for calibrating i?23- 

Vflchez & Esteban (1996) and Diaz & Perez- Montero (1999) have explored the use of the sum of 
the [S II] and [S III] forbidden-line strengths, S 23 = ([S II] AA6716, 6731 + [S III] AA9069, 9532)/H/3 
as an empirical abundance index, in analogy to R23. The main advantage of this index over R23 
is its monotonic behavior as a function of oxygen abundance, and its insensitivity to the global 
degree of ionization of the nebula, at least for giant extragalactic HII regions. The right panel of 
Fig. 7 shows the behavior of this index for the multiple position observations of 30 Dor, Carina, and 
Orion. In contrast to R23 (left panel), the values show much more scatter and nearly complete 
overlap between the three HII regions. This probably is due to a combination of observational 
error in the [S III] AA9069, 9532 measurements and large local variations in [S II] AA6716, 6731/H/J, 
which forms preferentially at the ionization interfaces and in shocks. The S23 ratio (and hence 
the inferred metal abundance) also varies systematically across the nebulae, especially in 30 Dor, 
where the data span the full range of radius in the HII region. This suggests that S23 should 
only be applied to global spectra of HII regions, which sample a representative fraction of the 
S ++ and S + emitting volumes. Despite the considerable scatter in 523 within each HII region, 
the corresponding dispersion in inferred metal abundances is only slightly larger than for R23, due 
to the relatively weak Z-dependence of S23 in this excitation range, if the calibration of Diaz & 
Perez-Montero (1999) is used. 

Figure 8 shows a similar test of the [O III]/[N II] empirical abundance parameter. We would 
expect this ionization-sensitive index to show a larger dispersion within individual HII regions, 
and this is borne out in Fig. 8. The index shows a dispersion of >1 dex in all three HII regions, 
and the inferred abundance values overlap between the objects. The horizontal lines indicate the 
corresponding oxygen abundances in terms of 12 + logO/H, using the calibration of Dutil & Roy 
(1999); the inferred abundances show dispersions of ±0.2 — 0.3 dex in a given HII region, about 2-3 
times larger than for R23. The relevance of these results for extragalactic applications is unclear, 
however. Comparisons of oxygen abundances derived from [O III]/[N II] and i?23 have tended to 
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show excellent agreement across a wide range of excitation (Ryder 1995, Roy & Walsh 1997), and 
indeed the former abundances show a lower scatter, which suggests that the [O III]/N II] method 
can be reliable for luminous HII regions where the spectrophotometric aperture encloses most of 
the emitting volume. However it is also possible that the apparent consistency of the [O III/N II] 
abundances is the result of a relatively narrow range of ionization parameters and other nebular 
properties within the samples that have been studied to date, and that this consistency could break 
down when applied to samples of HII regions with systematically different aperture sampling, 
ionization structure, or other nebular properties. 

4.3. Diagnostics of Stellar Temperature and IMF 

Emission-line spectra of H II regions and starbursts are being applied increasingly to constrain 
the ages and IMFs in the exciting star clusters (e.g., Stasihska & Leitherer 1996, Garcia- Vargas 
et al. 1997, BKG, Garnett 2000, and references therein). These applications rely on the ability 
to use nebular ionization models to characterize the shape and hardness of the ionizing continuum 
shortward of the Lyman break. In this section we use the sample of calibrating HII regions to 
directly assess several hardness parameters, ranging from robust indices such as the ionized He 
fraction to crude but commonly used indices such as the [O IIl]/H/3 ratio. 

4-3.1. He Recombination Lines 

The fractional ionization of He provides a robust measure of stellar temperature, following 
methods established originally by Zanstra (1927). Most of the helium in normal H II regions is either 
neutral or singly ionized, so the relevant diagnostics are the ratios of the He I and H I recombination 
lines. Figure 9 shows the line ratios He I A5876/H/3 and He I A6678/Ha plotted as functions of 
stellar temperature. Superimposed are the model relations, based on the stellar atmosphere models 
of BKG and assuming Case B recombination following Osterbrock (1989). Models are plotted for 
abundances of 1.0 Z Q (solid lines) and 0.2 Z (dashed lines), and for ionization parameter \ogU of 
—1.5 (upper) and —4.0 (lower). Given the relatively narrow range of U in our sample the ionization 
parameter dependence is unimportant, but the range of abundance in the sample has some effect 
on the observed line ratios (both via line blanketing in the ultraviolet and the He abundance) . The 
symbols differentiate between calibrating single-star regions, and larger Galactic and Magellanic 
Cloud HII regions, following Figure 4. The model lines represent the mean He ionization over the 
entire nebula, so for the H II regions with partial He ionization we only considered objects for which 
we had full or representative spatial coverage of the nebula. 

Figure 9 shows an excellent agreement between the observed levels of He ionization and the 
model predictions, especially for He I A6678/Ha where our data are the best and reddening errors 
are minimal. The most discrepant regions are S275, where our spatial coverage is very incomplete, 
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and S162 (NGC 7635), which is ionized by an 06.5 Illf star. The latter is the only HII region in our 
sample which is predominantly ionized by a giant star, but it is unclear whether the discrepancy 
for this object is coincidental or indicative of a general problem with the low-gravity stellar models. 
Otherwise the data show that the He I lines provide an excellent measure of stellar temperature. 
Unfortunately this method is only sensitive when nebular He is partially ionized (35000 K < 
%ff < 39000 K). Most extragalactic HII regions show nearly complete ionization of helium (%g 
> 40000 K), but the He I lines do provide important constraints on the ionizing spectrum in 
the innermost metal-rich disks of spirals, where many other hardness indices break down (BKG). 
Moreover, Fig. 9 provides independent confirmation of the reliability of the Vacca et al. (1996) 
effective temperature scale, at least for luminosity class V stars in the %g = 34000 — 39000 K 
range. 



4-3.2. Forbidden- Line Excitation Indices 

The simplest and crudest indicators of stellar temperature are based on single forbidden-line 
ionization indices, such as [O IIl]/H/3 (Chopinet & Lortet-Zuckermann 1976, Kaler 1978, Copetti, 
Pastoriza, & Dottori 1985, 1986), or analogous infrared indices such as [O III] 88 /im/H53a or 
[N III] 57 jum/H 53a (Puxley et al. 1989). Figure 10 shows the dependence of [O III] AA4959, 5007/H/3 
on stellar temperature, again with the H II regions subdivided as in Figure 4. This measure of the 
excitation shows a tight correlation with %g (for this sample of HII regions), confirming the good 
correlations seen previously by Chopinet & Lortet-Zuckermann (1976), Kaler (1978), and Hunter 
(1992). A roughly linear, increasing trend trend in [O III]/H/? with %g is expected, because the 
nebular volume containing double-ionized oxygen increases with the higher ionizing flux emitted by 
earlier spectral type stars (Stasihska 1978). The observed correlation has a remarkably linear form, 
extending to the WR stars with temperatures of nearly 70000 K. The adopted %]g scale for the 
WR stars (Esteban et al. 1993) is partly based on nebular modeling, however, so strictly speaking 
the empirical relation in Fig. 10 only applies for %g < 50000 K. 

Despite the tightness of the correlation in Figure 10, the [O IIl]/H/3 ratio must be applied with 
considerable caution to extragalactic HII regions and starbursts, because the ratio is sensitive to 
other variables such as ionization parameter and metal abundance. This is illustrated by the four 
model sequences in Fig. 10, plotted for oxygen abundances of 0.2 Z & and 1.0 Zq and log U = —2.0 
and —4.0 (with %g varying along each sequence). The best fitting models correspond to typical 
ionization parameters of log?7 ~ — 2.5 and Z/Zq ~ 0.5 — 1, consistent with the measured values 
from Fig. 5 and Table 1, respectively. However the same excitation vs %g relation may not apply 
in external galaxies, where both Z/Z & and logU may be shifted relative to the typical values in 
this sample. 

We have also investigated the behavior of two other ionization indices which have been sug- 
gested as hardness indicators, [Ar III] A7135/Ha and [Ne III] A3869/[0 II] A3727. Both are plotted 
against stellar temperature in Figure 11. The figure also shows the temperature dependence of 
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[O III] AA4959, 5007/H/3, but plotted this time with a logarithmic scale to provide a direct com- 
parison with the other two indices. The behavior of [Ar III] /Ha is qualitatively similar to that 
seen for [O IIl]/H/3, but with a larger scatter and shallower slope above %ff~ 38000 K. This is 
not surprising because the lower ionization potential for Ar + (27.6 eV vs 35.1 eV for + ) is close 
to the He + edge at 24.6 eV. Furthermore the excitation of [Ar III] is sensitive to the ionization 
parameter and abundance in the same way as [O III] , so it provides a poor substitute for the much 
more robust He + /H + index discussed earlier. 

The bottom panel of Fig. 11 shows the temperature dependence of the hybrid forbidden-line 
ratio [Ne III] A3869/[0 II] A3727. Ali et al. (1991) have shown that this ratio provides a valuable 
index for deriving ionization correction factors in nebular abundance measurements. Figure 11 
shows a strong monotonic dependence of [Ne IIl]/[0 II] on T^ff, which reflects the 41.1 eV ionization 
potential of Ne + , but the large scatter in the relation limits its usefulness as a stellar thermometer. 
The weakness of the [Ne III] line in metal-rich HII regions also limits its value for extragalactic 
applications. A similar plot of [Ne IIl]/H/3 shows a nearly identical scatter, so variations on [O II] 
excitation are not responsible for this dispersion. Our result is somewhat surprising in light of the 
relatively robust behavior of [Ne IIl]/H/3] in the recent analysis of Oey et al. (2000). Of all of 
the forbidden-line ratios that are accessible at visible wavelengths, the familiar [O IIl]/H/3 index 
appears to be the most useful stellar temperature indicator, though this sensitivity is accompanied 
by strong dependences on many other nebular properties. 



4-3.3. Composite Ionization Indices: rf 

The well-known sensitivity of excitation indices such as [O III] /H/3 to nebular ionization struc- 
ture and abundance led Mathis (1985) and Vilchez & Pagel (1988) to investigate the use of compos- 
ite indices based on the ionization of more than one element, in an effort to find a more robust UV 
hardness indicator. The advent of CCD spectrographs has made the [S III] AA9069, 9532 doublet 
accessible, and with it the possibility of using the + /0 ++ and S + / S ++ to contrain %g indepen- 
dently of metal abundance and ionization parameter. A convenient hardness index is the parameter 
rf introduced by Vflchez & Pagel (1988): 

_ [O J/JAA3726, 3729/[Q III] AA4959, 5007 
V ' ~ [S II] AA6716, 6731/ [S III] AA9069, 9532 ' ^ ' 

The robustness of rf can be tested directly by plotting the behavior of [O II] /[O III] vs [S II] /[S III] 
for 30 Doradus, Orion, and Carina, where we have spectra at multiple positions. This is shown 
in Figure 12a, along with a series of nebular models from BKG. Constant %^ lines are drawn at 
Z = 0.5 Z Q for % = 50000, 45000, 40000, and 35000 K (left to right). Along each line the ionization 
parameter varies from logU = —1.5 (bottom) to logU = —4.0 (top). The dashed lines show the 
effects of changing metal abundance between Z = 0.2 Z Q (left) and 1.0 Z Q (right, indicated by small 
dots), for the same range of stellar temperatures and a fixed ionization parameter logU = —2.5. 
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For further details on the models see BKG. 

For 40000 K < %g < 50000 K the lines of constant stellar temperature have a slope of near unity 
in this diagram (i.e. constant rj'), demonstrating the robustness of rj against changes in ionization 
parameter, at least in the models. More impressive yet is the distribution of the observed line ratios 
at the various positions in 30 Doradus (open stars) and Orion (open circles), which fall on lines of 
nearly constant r/, over ranges of nearly an order of magnitude in [O II] /[O III] or [S II/ [S III]. The 
observations of Carina (triangles) show a large scatter in rj', and we suspect that this dispersion is 
real, reflecting the influence of ionization from stars with different temperatures across this complex 
region. However Oey et al. (2000) have also found that rf shows a larger point-to-point variation 
in their sample of 4 LMC shell HII regions, and this may suggest limits to the range of spatial 
scales and nebular properties where rf provides a reliable measure of ionizing stellar temperature. 

Although the spectra of 30 Dor and Orion are consistent with a nearly constant value of rj' (and 
hence %g), Fig. 12a shows a systematic offset between the effective temperatures implied by the 
models (~39000 K and 37000-38000 K respectively), and those estimated from the stellar contents 
of the regions (~48500 K and 40000 K respectively). Figure 12b shows the same diagram for the 
other HII regions in the sample (a nebular model for %g = 30000 K is added), and a similar offset 
in the stellar temperature scales is apparent, if one compares with the actual %g values in Table 1. 
The HII regions ionized by cooler stars {%g < 37000 K) seem to show better agreement with the 
models. Nevertheless the general consistency between the nebular and spectroscopic temperatures 
is encouraging, confirming and extending earlier results by Mathis &; Rosa (1991). 

Figure 12 impressively demonstrates the robustness of rj' against systematic variations in U , 
but the models plotted in the figure suggest a significant dependence of rj' on metal abundance. 
The dashed lines show the effects of changing Z from 0.2 Z (left) to 1.0 Zq (right). The in- 
dependence is relatively weak for hotter stars (%g > 40000 K) and low abundances (Z < 0.5 Z ), 
but becomes large for cooler stars, where the effects of abundance and stellar temperature become 
virtually degenerate, and rj' itself becomes nearly degenerate with [O II/O III] and [O IIl]/H/3. 
Inspection of the models suggests that two physical effects contribute to this abundance sensitivity, 
increased stellar metal line blanketing in the ultraviolet, which softens the ionizing spectrum at 
a given effective temperature, and the increasing dominance of nebular cooling (especially ++ ) 
over ionization structure in determining the relative forbidden-line strengths at high abundance 
(Shields & Searle 1978, Oey & Kennicutt 1993, Stasinska & Leitherer 1996). Thus it is important 
to confirm the ?/-derived stellar temperatures for low-excitation nebulae with those derived from 
more robust hardness indices, such as the He recombination lines. 

The degeneracy between abundance and stellar temperature in model spectra has hampered 
previous model-based studies of the metallicity dependence of the stellar temperatures and IMFs 
of extragalactic HII regions (e.g., BKG). However the data presented in this paper allow to test 
the correlation between rj' and %g empirically, independently of the nebular models. The results 
are shown in Figure 13; the points show the measured values of log rj' and %g , subdivided as in 
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Fig. 4, while the three lines show the nebular model relations for 0.2, 0.5, and 1.0 Z Q (log U = —2.5 
in all three cases). The HII regions span a range of roughly 2 dex in log rj (S237, with no 
detected [O III] emission, is not plotted), and the index shows a monotonic dependence on stellar 
temperature, albeit with a considerable scatter. The general shape of the log 77/ vs %g relation is 
consistent with that predicted by the nebular models; the hatched line shows a quadratic fit to the 
observed HII regions (excluding the WR nebulae, where the temperature calibration is uncertain), 
and its general shape is in good agreement with the model relations for the relevant abundance 
range of ~0. 4-1.0 Z . 

In contrast to the behavior of our nebular models, our HII region data show no evidence for 
a strong Z-dependence of the 7/ vs %jj relation. In particular, there is no evidence in Fig. 13 for a 
systematic shift between the LMC HII regions and the more metal-rich (on average) Galactic HII 
regions, or likewise for any difference between the most metal-rich and metal-poor Galactic HII 
regions. Within the considerable scatter the trends seen in Fig. 13 are independent of abundance. 
The singular exception of the SMC HII region N66, which is the most metal-poor region in the 
sample, but the %g value of this region almost certainly is underestimated, due to the presence of 
the luminous WR star HD 5980 (see the appendix). 

We cannot explain why the observed behavior of rf appears to be more robust than the models 
would predict, but we can speculate on the possible reasons. As pointed out earlier, our empirical 
estimates of %g for the HII regions are based on stellar spectral classifications, using a single, 
metallicity-independent spectral type vs %g calibration (Vacca et al. 1996). We believe that this 
procedure is justified, because the stellar types are predominantly assigned on the basis of H, He , 
and He + absorption line strengths, which ought be robust against metallicity effects. However 
the stellar atmosphere grids which are used in the nebular models do show considerable changes 
in ultraviolet line blanketing with increasing metal abundance, and this may alter the ionization 
temperature at a given fixed effective temperature, in a way that would not be reflected in the 
assigned spectral types of the stars. Stated in another way, although we have calibrated 7/ in terms 
of stellar effective temperatures in Figure 13, the ionization temperatures that are derived from 
7/ may not necessarily coincide with %g for abundances that differ substantially from the local 
Galactic value. We return to this point in §5. 

In summary, our data suggest that 7/ offers the prospect of estimating the characteristic stellar 
temperatures of the ionizing associations in HII regions, over a relatively large temperature range 
(33000-fT < T e ff < 50000-fT), and on a calibration that is directly tied to the stellar spectral 
classification system, independent of the nebular model chain. However two important limitations 
of 77' must be recognized before the results in Fig. 13 are applied indiscriminantly to HII regions 
in external galaxies. First, the scatter in the log 77' vs %g relation is considerable, amounting to 
an rms uncertainty of ±3000 K for any individual region. Second, photoionization models suggest 
that abundance effects may introduce systematic shifts of a comparable magnitude in the derived 
effective temperatures, though the observations of the calibrating HII regions do not show evidence 
for shifts of this magnitude. 
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5. Application to Extragalactic HII Regions 

In this section we briefly apply the empirical relations between nebular hardness properties 
and stellar temperature from the previous section to the sample of extragalactic H II regions stud- 
ied by BKG. The BKG data are well suited to our purposes because they span a large range of 
metal abundance and spiral galaxy type, and they are based on spectra covering the 3600-9700 A 
range, so we can apply all of the hardness indices discussed earlier, including rj' . Our approach is 
complementary to the model-based studies of Vflchez & Pagel (1988), Stasihska & Leitherer (1996), 
and BKG. 

We begin by using the empirical fit to the log rf vs %g relation in Fig. 13 (hatched line) to 
derive stellar temperatures for the BKG extragalactic HII regions with reliable measurements of 
rj'. Fitting a quadratic function to the calibrating Galactic HII regions (solid squares) yields: 

T eff = 50819 - 16485a; + 3778a; 2 (3) 

where x = log rf. Although this relation is based solely on the smaller calibrating HII regions, 
where the %g values are the most reliable, using the entire sample does not change the resulting 
relation significantly. Equation (3) is given primarily so readers can reproduce and check the results 
described below; in view of the caveats discussed in the previous section we strongly advise against 
applying this relation generally to extragalactic regions, least of all to luminous starbursts, where 
the physical conditions may differ substantially from those of the calibrating HII regions. 

With these caveats firmly in mind, the top panel of Figure 14 shows the characteristic stellar 
temperatures T* and oxygen abundance for the BKG sample. We have deliberately avoided referring 
to these as effective temperatures, to emphasize that these are ionization temperatures which may 
not coincide precisely with %g outside the calibrated abundance range. The abundances have 
been derived using the i?23 empirical method following the calibration of Edmunds & Pagel (1984), 
and should be reliable to within ±0.2 dex for Z < Z , but less accurate at higher abundances. 
However this accuracy is sufficient for testing for general trends in T* with metallicity. The values 
of T* inferred from equation (3) show a large scatter, consistent with the large scatter observed 
in the calibrating sample (Fig. 13), combined with any real scatter in T* at fixed abundance, due 
to variations in age and other properties of the ionizing star clusters. The large dispersion in Fig. 
14a underscores the dangers of applying rj' or any other forbidden-line hardness parameter to an 
individual H II region. Nevertheless the general behavior of T* should be meaningful. 

The top panel of Fig. 14 shows a clear trend in stellar temperature, with the typical %g 
decreasing from ~55000 K at 0.2 Z Q (12 + logO/H = 8.2)to ~40000 K above solar abundance 
(12 + log O/H = 8.9). This trend is in excellent agreement with the model-based results of BKG, 
except for a temperature offset of about 5000 K; BKG esitmated a %g range of 50000 K to 35000 
K over the same abundance range. In both cases note that most of the %g change is restricted to 
abundances below solar. 

The measured He I line strengths in the BKG spectra provide another empirical constraint on 
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the stellar temperatures in metal-rich HII regions. Unfortunately the He lines are only useful for 
regions ionized by stars with %g < 38000 K, when one allows for observational uncertainties in 
the He line fluxes and the dependences on abundance and ionization parameter. Nevertheless we 
can at least test the consistency of the rf and He + /He° temperature scales for the most metal-rich 
regions (12 + log O/H > 9.0). The crosses in Fig. 14 show the values of %g that are implied by the 
observed He IA5876/H/3 ratios and the models shown in Fig. 9 (the He IA6678 line is usually too 
faint to provide reliable constraints on %g), with with arrows indicating objects for which the He 
recombination lines only set lower limits on %g. The He I data confirm the general trend toward 
lower T* for Z > Zq, but with T* ~ 38000 — 39000 K at the highest abundances, about 2000 K 
cooler than given by rj'. This comparison provides a good illustration of how two different hardness 
indices can provide complementary constraints on the stellar populations in the HII regions. 

As an additional consistency check on these results, the bottom panel of Fig. 14 shows the same 
relation, but in this case with the T* values derived from the average of the rf vs %g calibration 
(Fig. 13 and eq. [3]), and the relation between [O IIl]/H/3 and %g from Fig. 10). The latter is 
known to be sensitive to ionization parameter and to metal abundance itself, but surprisingly it 
yields values of T* that are close to those derived from rf . This comparison yields a somewhat 
lower T* values at high abundance, which are in better agreement with the values derived from the 
He I recombination lines. However because of the susceptibility of [O IIl]/H/3 to dependences on 
parameters other than %g we do not apply this diagnostic further. We suspect that the general 
consistency of the results here are due to the fact that mean abundances and ionization parameters 
of the calibrating H II region samples are typical of the extragalactic H II regions in the BKG sample, 
but this does not imply that the abundance dependence of T* implied by the [O III] excitation is 
necessarily robust by any means. 

In summary the empirically-based stellar temperature indices show evidence for a systematic 
decrease in mean stellar temperature with increasing metal abundance, which appears to quali- 
tatively confirm previous results based on nebular model fitting (Vilchez & Pagel 1988, BKG). 
Several physical effects can contribute to this softening of the ionizing continuum with increasing 
metallicity: increased line blanketing at constant %g, a shift in %g at constant stellar mass, and/or 
a change in the upper mass limit to the IMF. The models of Stasinska & Leitherer (1996) suggest 
that changes in atmospheric blanketing and stellar evolution can account for most of the changes, 
and our results do not contradict this conclusion. We refer the reader to that paper and to BKG 
for a more complete discussion of the implications of these changes in %g . 

6. CONCLUSIONS 

In conclusion, we have presented a set of empirical diagnostic line ratios, calibrated against 
the effective temperature of the exciting stars. Our results provide a stronger empirical foundation 
for many of the most widely used nebular diagnostics (e.g., R23, rj'), and they show much better 
agreement between observations and models than often was seen in earlier generations of stellar 
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atmosphere and nebular photoionization models (cf., Mathis 1985, Vilchez k, Pagel 1988). The 
empirical relations presented here should be useful in many investigations of the stellar populations 
of extragalactic HII regions, where very limited information is available on the exciting stars, and 
where an estimate of the nebular-based effective temperature is needed, even with a restricted 
wavelength coverage. 

We close by reemphasizing the exploratory nature of this study, and pointing out the most 
urgently needed improvements before one can reliably extract quantitative constraints on the stellar 
contents and IMFs of star forming regions from nebular spectra. On the empirical side, this 
approach can be much improved if accurate (T e -based) metal abundances are obtained for all 
of the calibrating HII regions (or their exciting stars), and the effects of effective temperature, 
metal abundance, and nebular geometry are accounted for explicitly in the analysis. Oey et al. 
(2000) have demonstrated the value of obtaining high signal/noise point-by-point measurements 
of individual H II regions, and applying this approach to a larger sample of objects would provide 
a major improvement over our exploratory analysis. On the theoretical side, more self-consistent 
stellar models for the relevant range of abundances, temperatures, and gravities are the most crucial 
inputs. This would also allow for a consistent treatment of stellar and nebular abundances in the 
photoionization calculations. We hope that these initial results may stimulate further progress in 
these areas. 

It is a pleasure to thank A. Pauldrach, R.-P. Kudritzki, and T. Hoffmann for sharing their 
unpublished stellar atmosphere models, and S. Oey for a preliminary version of her paper and 
many discussions about this work. This research was supported in part by the National Science 
Foundation through grants AST-9421145 and AST-9900789. 

APPENDIX 

EXPLANATORY NOTES ON INDIVIDUAL H II REGIONS 

S257: The spectral type adopted is the one predicted by HM90 on the basis of the observed total 
nebular ionizing flux. 

S275 (Rosette): The ionizing cluster is NGC 2244; earliest spectral type is 04 V((f)); a total of 9 
stars earlier than BO. 5 has been included in the computation of %g. 

RCW 5 (NGC 2359): A Wolf-Rayet (WN 4) ring nebula. Esteban etal. (1993) derive %=67±11 
xlO A K. Given the uncertainties in the treatment of WR star atmospheres, this object, together 
with RCW 48, is included to observe the trends in the nebular diagnostics at very high %g values. 

RCW 16: Walborn (1982) gives two early-type stars in NGC 2467, HD 64568 (03 V((f*))) and HD 
64315 (06 Vn). 

RCW 48 (NGC 3199): The second Wolf-Rayet (WN 5) ring nebula in our sample, with %=57±14 
xW A K (Esteban etal. 1993). 
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RCW 53 (Carina): A total of 36 early-type stars in Trl4 and Trl6 has been considered, excluding 
r\ Carinae and one WR star (HD 93162). 

RCW 57 (NGC 3603): Drissen et al. (1995) identify 11 O-type stars and 3 WN6 stars which account 
for at least 80% of the ionization of the nebula. The temperature derived here does not include the 
WN6 stars, which according to Esteban et al. (1993) have % s < 42000 K. 

S49 (M16): 20 early-type stars 

S45 (M17): 12 early-type stars classified from optical and if-band spectroscopy by Hanson, Howarth 
& Conti (1997). 

S100: The spectral type adopted is the one predicted by HM90 on the basis of the observed total 
nebular ionizing flux. 

S158 (NGC 7538): Ionizing star is CGO 654, with spectral type 07V, from Crampton, Georgelin, 
& Georgelin (1978). 

S162 (NGC 7635): Ionizing star is BD+60 2522, with spectral type 06.5 Illf, from Conti & Alschuler 
(1971). Classified as 06.5 Illef by Conti & Leep (1974). 

N11B: We have assumed that the stellar association LH10 is responsible for the ionization (26 
ionizing stars). 

N44: The slit was centered on N44B, ionized by LH 47/48. 

N51D: We have included the 40 early-type stars classified by Oey k, Smedley (1998) in LH 51 and 
LH 54. 

N144: Garmany, Massey & Parker (1994) list 32 ionizing stars in LH 58. The 3 WR stars were not 
included. 

30 Dor: We have included the 64 early- type stars in R136 listed by Massey & Hunter (1998), 
excluding the WR stars. 

N70: The ionizing cluster is LH 114, studied by Oey (1996). 

N180: The ionizing cluster has been assumed to be LH 117, studied by Massey et al. (1989a). 

N66 (SMC): The temperature we derived is based on the stellar census of Massey et al. (1989). 
However the calculation does not include the peculiar WN3+OB binary HD 5980, which according 
to Massey et al. could account for up to half of the ionization of the nebula. Consequently the 
temperature listed in Table 1 is strictly a lower limit. 
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FIGURE LEGENDS 

Fig. 1. — Distribution of oxygen abundances and ionization- weighted stellar temperatures for the 
HII region sample. Open squares denote nebulae with direct abundance determinations based 
on electron temperature measurements. Filled triangles denote HII regions with empirical (R23) 
abundance estimates. 

Fig. 2. — Charts showing the locations of the long-slit spectrophotometric measurements. In all 
cases north is at the top and east is to the left. The box in M8 (Fig. 2d) indicates the region 
covered by the drift scan. 

Fig. 3. — Diagnostic diagrams showing the differences between the fixed pointing spectra (solid 
points) and the drift scans (open points) for the northern sample of Galactic H II regions. 

Fig. 4. — Comparison of spectral properties of the Galactic calibrating HII regions (solid squares), 
large Galactic HII regions (open squares), LMC and SMC HII regions (open circles), and WR neb- 
ulae (open stars). The small triangles show for comparison the spectral sequence for extragalactic 
H II regions in a large sample of galaxies from BKG. 

Fig. 5. — Relation between [O IIl]/H/3 and [N II] /Ha for the HII regions in this sample, subdivided 
by stellar effective temperature (left) and by nebular oxygen abundance (right). 

Fig. 6.— The [S II]AA6716,6731/[S III]AA9069,9532 vs % diagram for the HII regions. The symbol 
coding is the same as in Fig. 4. The continuous lines show nebular models at 0.2 Z@, for a variety 
of ionization parameter values log U (as labelled). The dotted lines refer to solar abundance models 
at log U = -2.0, -2.5 and -3.0. 

Fig. 7. — (Left): Relation between [O IIl]/H/3 and [O Il]/H/3 for subregions observed in 30 Doradus 
(open stars), Carina (solid triangles), and Orion (open circles). The dashed lines show constant 
values of the empirical abundance parameter R23. Note the relative constancy of R23 across the 30 
Dor and Carina, despite the large local variations in excitation. (Right): Corresponding relation 
between [S IIl]/H/3 and [S Il]/H/3 line ratios. Diagonal lines show constant values of the abundance 
index 523. 

Fig. 8. — Relation between [O IIl]/H/3 and the extragalactic HII region empirical abundance index 
[O III]/[N II], for subregions observed in 30 Doradus (open stars), Carina (solid triangles), and 
Orion (open circles). The dashed horizontal lines show values of constant oxygen abundance, using 
the calibration of Dutil k Roy (1999). 

Fig. 9. — Dependence of nebular He I recombination line strengths (relative to Balmer lines) as a 
function of stellar effective temperature, as measured by He IA5876 (top) and He IA6678 (bottom). 
Both ratios are approximately proportional to the fraction ionization of He in the nebula. Symbols 
are coded as in Fig. 4. The lines show the relations expected from stellar atmosphere and nebular 
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photoionization models for different abundances, as described in the text. 

Fig. 10. — Correlation between excitation index [O III]/H/? and stellar effective temperature, with 
symbols coded as in Fig. 4. The four sets of lines show photoionization model sequences for oxygen 
abundances of 0.2 and 1.0 Z & and ionization parameters \ogU = —2.0 and —4.0. 

Fig. 11. — Dependences of three nebular excitation indices on stellar effective temperatures. The 
symbols are coded as in Figs. 4 and 10. 

Fig. 12. — Relationship between ionization-sensitive forbidden-line ratios [S II] /[S III] and 
[O II] /[O III], (a) Spatially resolved observations of the 30 Doradus nebula (open circles), the 
Orion nebula (solid squares), and the Carina nebula (open triangles). The solid lines show pho- 
toionization model sequences for stellar temperatures of 50000 K, 45000 K, 40000 K, and 35000 K 
(left to right), (b) The same line ratios for the full sample of HII regions, with symbols coded as 
in Fig. 4. Solid lines show the same model sequences, with the addition of a 30000 K model (far 
right). Dashed lines show the effects of increasing the metal abundance, as discussed in the text. 

Fig. 13. — Dependence of the ionization hardness parameter rf (eq. [2]) and stellar effective for 
the HII regions in this sample, with symbols coded as in Fig. 4. The lines show photoionization 
model sequences for oxygen abundances of 0.2, 0.5, and 1.0 Z & . The hatched line shows a quadratic 
empirical fit to the data. 

Fig. 14. — The top panel shows the relationship between mean stellar temperature inferred from 
the rj' parameter (eqs. [2, 3]) and nebular oxygen abundances, for the BKG sample of extragalactic 
HII regions. The middle panel shows the stellar temperatures derived from an average of r/ and 
[O IIl]/H/3, as described in the text. Crosses in both panels show values of T* derived from He I 
A5876/H/3. Arrows indicate lower limits to T*, corresponding to full He ionization. 



-27- 



Table 1. H II region sample 
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References. — Sources of spectral types: 1. Hunter & Massey 1990; 2. Conti & Alschuler 1971; 3. Massey etal 1995; 4. 
Lahulla 1987; 5. Esteban etal 1993; 6. Walborn 1982; 7. Heydari-Malayeri 1988; 8. Georgelin etal 1973; 9. Massey & Johnson 
1993; 10. Drissen etal 1995; 11. Walborn 1987; 12. Lada etal 1976; 13. Hillenbrand etal 1993; 14. Hanson etal 1997; 15. Parker 
etal 1992; 16. Oey & Massey 1995; 17. Oey & Smedley 1998; 18. Garmany etal 1994; 19. Massey & Hunter 1998; 20. Oey 
1996; 21. Massey etal 1989a; 22. Massey etal 1989b. 

Sources of abundances: a. Pcimbert etal 1993; b. This paper; c. Mclnick etal 1989; d. Shaver etal 1983; e. Talent & Dufour 
1979. 



-28- 



Table 2. Reddening-corrected line fluxes 
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